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Abstract 

If j4 is a set of nonnegative integers containing 0, then there is a unique nonempty set B of 
nonnegative integers such that every positive integer can be written in the form a + b, where 
a G A and 6 £ _B, in an even number of ways. We compute the natural density of B for several 
specific sets A, including the Prouhet-Thue-Morse sequence, {0} U {2": n € N}, and random 
sets, and we also study the distribution of densities of B for finite sets A. This problem is 
motivated by Euler's observation that if A is the set of n that have an odd number of partitions, 
then B is the set of pentagonal numbers {n(3n+ 1)/2 : n € Z}. We also elaborate the connection 
between this problem and the theory of de Bruijn sequences and linear shift registers. 



1 Introduction 

There is a unique set B of nonnegative integers with the property that each positive integer can be 
written in the form + 6 (s G N := {0, 1, 2, . . . }, 6 G i?) in an even number of ways. Specifically, 

B = {0, 1, 2, 3, 5, 7, 8, 9, 13, 17, 18, 23, 27, 29, 31, 32, 35, . . .}. 

Are the even numbers in B exactly those of the form 2fc^? Does B have positive density? 

Before addressing these two questions, we restate and motivate the problem in greater generality. 
Given any sets A, i? C N, the asymmetric additive representation function is defined by 

R{n) #{(a, 5) : n = a + 6, a e A, 6 e B}; 
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Figure 1: The points {n, |Podd H [0,n]| - jPoven n [0,n]| - l) = (n,2|Podd n [0,n]| - 



n 



equivalently, we could define R by noting that 



beB 



We are interested in the situation where i?(0) = 1 and R{n) = (mod 2) for n > 0, i.e., the 
situation where R{'n)<f^ = 1 in the ring of power series F2[[q]]. In this case, we say that A and 
B are reciprocals, and we write A = B and B = A. The general problem of this paper is to find 
the reciprocals of several special sets A, and to draw some conclusions about "typical" properties of 
reciprocals. We are particularly concerned with the relative density. 



5{A, n) 



|An [0,r 



and the density 5{A) := liuin^ac ^{A, n) (when the limit exists). 

We began studying this problem after reading two articles by Berndt, Yee, and Zaharescu [SJE], 
where bounds on the density of the set Podd := {n e N: p{n) = 1 (mod 2)}, with p{n) being the 
ordinary partition function**, are proved. The starting point for their work is Euler's pentagonal 
number theorem |12l Theorem 10.9], and in particular that the reciprocal of Podd is the set {n(3n + 
l)/2: n g Z} of pentagonal numbers. Since the known bounds (see and JT] for the currently-best 
results) on the thickness of Podd are so strikingly far from what is believed to be true, we felt that 
it would be beneficial to study the "reciprocal" notion in a more general setting. 

^odd is pictured in Figure ^ where not only does it appear to have density 1/2, but the walk 
defined by w{n) := 2 |Podd H [0,n]| — n visually appears to be a simple random walk. See |18j for 

''p(n) is the number of ways to write n as a sum of nonincreasing positive integers. For example, 4 = 3 + 1 = 
2-f2 = 2-|-l + l = l + l + l + l, so p(4) = 5. 
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a report of more elaborate statistical tests on the set Podd- We note that while p{n) appears to be 
uniformly distributed modulo 2 and 3, it has been known since the time of Ramanujan to not be 
uniformly distributed modulo 5, 7, or 11. 

In contrast to that of the pentagonal numbers, the density of the reciprocal of the squares 
appears to drop off steadily to 0. Set S :— {in? : n e N}, with reciprocal S. The relative density of S is 
pictured in FigureEl In Section f6.il we prove that the even numbers in S are precisely {2n^ : n e N}, 
and we characterize the n G S with n = 1 (mod 4) as those n whose prime factorization has a 
particular shape. Those n e 5 with ri = 3 (mod 4) are characterized in terms of the number of 
representations of n by certain quadratic forms. 

Generalizing the squares and pentagonal numbers, we treat 



e(ci,C2) := <^ Cin + C2 



n{n — 1) 



n e Z 



for general c\ and ci in Section El A few interesting special cases are the binomial coefficients 
0(0, 1) — {(2) : n G N}, the squares 8(1, 2), and the pentagonal numbers 8(1, 3). 

Conjecture 1.1. The reciprocal of the set 8(ci,C2), where < 2ci < C2 and gcd(ci,C2) = 1, has 
density i!/c2 = 2 (mod 4), and otherwise has density 1/2. More precisely, if C2 = 2 (mod 4), then 



lim 

n — *oo 



e(ci,c2)n[o,n] 



nj logn 

for some positive constant C depending only on C2. //c2 ^ 2 (mod 4), then 



lim sup 



8(ci,C2) n [0,n] 



loglog(7i)/2 



= 1. 



Numerically, it seems that the constant C is 2 if C2 = 2 or 6, and C = 4 if C2 = 10. We lack 
sufhcient data to guess the other values. The authors believe that the non-effective C2 = 2 (mod 4) 
case might be provable by showing that the generating function of 8(ci,C2) is congruent modulo 
2 to an integer-weight modular form, which has almost all of its Fourier coefficients even. This is 
outside the scope of this paper, and we leave it as an area for further study. 

The C2 ^ 2 (mod 4) case is motivated by the celebrated law of the iterated logarithm. Let 
X\,Xi^ ... be independent random variables taking the values and 1 with probability 1/2. The 
law of the iterated logarithm states that 



lim sup 



v/nloglog(n)/2 



1 



with probability 1. What we actually would like to conjecture is that reciprocal of 8(ci,C2), with 
C2 ^ 2 (mod 4), is statistically indistinguishable from a truly random set with density 1/2. The 
phrase "statistically indistinguishable" is too vague, however, so in Conjecture II . II we have settled 
for this one specific statistic. 

The natural expectation is that, barring some cosmic coincidence or obvious structure, the 
reciprocal of a set should have density 1/2. This is affirmed by the case of a random set, which we 
handle in detail in Section |3| let X\^Xi-, ... be independent random variables taking the values 
and 1, with probabilities bounded away from and 1, and set F :— {0}U{n: X„ = 1}. Theorem l3.1l 
states that the reciprocal of F has density 1/2 with probability 1. This makes the sets whose 
reciprocals do not have density 1/2 the interesting ones. 
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Our purpose is to identify relevant properties of those sets whose reciprocals have density 
different from 1/2. Specifically, in addition to random sets and 6(ci, C2), we consider finite sets, the 
set of powers of two, and the set of Prouhet-Thue-Morse numbers. 

• Finite sets: the reciprocal has a rational density, and appears to typically have density slightly 
below 1/2. We identify through algebraic properties two infinite classes of polynomials, one 
whose reciprocals have density strictly larger than 1/2, and one whose reciprocals have density 
at most 1/2. 

• Powers of 2: the reciprocal of the thin set {0} U {2" : n e N} is the thin set {2" - 1 : 71 e N}. 
In particular, we describe the reciprocal of {0} U {2""' : rt g N} for every m e N. 

• Prouhct-Thue-Morse numbers'^: the reciprocal of 

T := {n E N: the binary expansion of n contains an even number of "l"s}. 

has density 1/3. Specifically, we prove that fc e T if and only ii k — or (fc± l)/4 is an integer 
whose binary expansion ends in an even number of zeros. 

The strongest conjecture that is consistent with our theorems, our experiments, and Coniecture ll.il 
is Conjecture II. 21 

Conjecture 1.2. // a set contains 0, is not periodic, and is uniformly distributed modulo every 
power of 2, then its reciprocal has positive density. 

We now include a section-by-section agenda for the remainder of the paper. 

Section ^ Motivate and contextualize reciprocals of sets. 

Section [21 Introduce notation and derive general expressions for reciprocals. 

Section |3l Consider reciprocals of random sets with positive density. 

Section 131 Consider reciprocals of finite sets. 

Section 151 Consider the reciprocal of the powers of 2, and similar sets. 
Sectional Consider the reciprocal of 0(ci,C2), particularly the squares. 
Section [3 Consider the Prouhet-Thue-Morse sequence. 
bT = {0, 3, 5, 6, 9, 10, 12, 15, 17, 18, 20, ... } 
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2 Notation and General Formulas 

Throughout this paper, we let 

Hq) = fo + fiq + f2q^ + --- and Hd) = fo + hq + hq^ + ■ ■ ■ (1) 
be elements of F2[[(i']] that satisfy the equation 

Hq)Hq) = 1- (2) 

In particular, fn — fo — 1. We define the integer sets F := {n>0: /„ = 1} and F :— {n > 0: fn — 

Note that ||2Jl implies (for all fc > 1) that J-{q^)J-{q^) — 1 also. This corresponds to noting that 
multiplying everything in by fc has the effect of multiplying everything in F by k. With this in 
mind, we sometimes make the convenient assumption that gcdi^ = 1. 

Our next lemma is a fundamental identity in F2[[q]], and has a number of remarkable conse- 
quences. We use it frequently throughout this paper. 

Lemma 2.1. The reciprocal of J-{q) is J-{q)J-{q'^)J-{q^)J-{q^) ■ ■ ■ . That is, 

oo 

i^T{q) ■ YlHq'^). (3) 

k=Q 

Proof. First, notice that both sides of this equation have constant term equal to 1. Also notice that 
for any fixed n > 0, only finitely many terms of the infinite product affect the coefficient of g". Thus, 
the coefficient of on the right hand side of Q is also the coefficient of q" in 

[log 2 "J 

Hq)- n 

fe=o 

By the so-called children's binomial theorem'^, T{q)J^{q) = ^{q'^)- Multiplying by !F{q^), we see 
that T{q)!F{q)!F{q^) = T{q^)T{q^) = J-{q^), and continuing we get 

Hq)Hq)Hq')Hq') ■ ■ ■ Hq'''°'' ) = Hi''"' • (4) 

Now notice that since < n < 2L'°S2 "J+i, the coefficient of q" on the right hand side of (@J is 0, 
and our result follows. ■ 

We now give a list of recurrences for /„, discuss the usefulness of each, and prove them. 
Lemma 2.2. If !F{q)J-{q) — 1, then fo = l and for n > 0, 

n 

ii- fn — ^ if o,nd only if ff^i {xq, xi, . . .) : Xi ^ F,n = y Xil^ \ is odd. 



Y^xA 

i>0 ^ 



Hi. fn = fxifx2 • • • fxi, where the summation extends over all tuples x = {xi,--- ,xe) with 

X 

n — X]i=i '^'^^ each Xi > (£ is allowed to vary); 



:(a + fe)2 = a2 + 62 (mod 2) 
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/n = X! fn-2ift + G(/i, /2, . . . , /l„/2J ), for some function G. 

0<i<n/4 

Lemma \2.2f i] is valuable because of its simplicity. For instance, it is immediately apparent 
from this recurrence relation that T is uniquely defined and always exists (provided fo = 1). 

In several of the examples we consider, the set F has some special properties modulo a power 
of 2. Lemma ij) facilitates our exploitation of these special properties. 

Lemma \2.2i iii ) is useful because of its symmetry, and because its right-hand side does not 
expressly reference the / sequence. As a specific example, let r{n) be the number of ways to write 
n as a sum of positive pentagonal numbers (counting order). Then, by Lemma l2.2r zM^. p{n) = r(n) 
(mod 2). We also use Lemma r2.2f wi). for example, to prove Lemma|^2Iw). 

If one lets the fi be independent random variables, then the expression in Lemma |^2Jw) 
contains a summation of weakly dependent random variables, and a deterministic function of 
/i, . . . , fnj2- This allows us to say something explicit about the distribution of the resulting random 
variable /„ (see Theorem 

Another remarkable aspect of Lemma \2.2i iv) is that /„ does not depend in any way on 
/„_i, fn~3, ■ • • , /n-c, where c is the largest odd number strictly less than n/2. For example, 

111 = fll + /9/1 + /7/2 + /7/1 + /5/3 + fifs + /3/2/1 + /3/2 + /2/1 + fl 

does not depend on /lo, /s, or /g. 

Proof. Comparing the coefficients of on the left- and right-hand sides of equation Q yields 
Sj=o fjfn-j = 0- Lemma l^^ z) is this expression rearranged, using the fact that /o = 1. 

Similarly, Lemma I^Jii) equates the coefficients of 9" on the left- and right-hand sides of 
equation (j^J, with the right-hand side interpreted as a product in Z. 

One can prove Lemma \2.2\ by induction, using Lemma \2.2\ i'\ to complete the induction 
step. Alternatively, one may simply compare the coefficients of on the left- and right-hand sides 
of 



which is valid because we are working over F2. 

Recall Kummer's result that the multinomial coefficient 



mi + ■ ■ ■ + m.k \ [mi + • • ■ + m^)! 



TOi,m2,...,mfe 



is relatively prime to a prime p if and only if toi, . . . , rrik can be added in base p without carrying |10| . 
We are working with p = 2, so our condition is: (™^^' ^'°) is odd if and only if no two of the binary 
expansions of mi, . . . ,mfe have a "1" in the same position. We call such a list of positive integers 
mi, . . . , m^, non- overlapping. 

Let 7r(n) be the set of partitions of n whose distinct parts xi,...^Xk have non-overlapping 
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multiplicities rni, . . . ,mk. Continuing from l2.2T zM). we have 

fn ^ ^ fxifx2 ' ' ' fx^ 



xi-\ yxi—n 

Xi>0 



^ ini -r • • ■ -r ink rm2 fm^. 



miaiH ^m^a^—n 

ai>->afc>0 
mi>0 

_ \ ^ p'mi rm2 _ , , j-rrifc 
/ J J a\ J a2 *' dfc 

7r(n) 

_ \ ^ rrni rm2 , , , mk i \ ^ fmi j-m2 _ , , j-mt 
/ y J ai J a2 •> ak ' / j J ai J a2 •' 

7r(n) 7r(n) 
ai>n/2 ai<?i/2 

If fli > n/2, then it must have multiplicity mi = 1, and if mi, . . . , nik are non-overlapping then the 
other m, are even: 



j?mi rm2 , , , j-nifc _ -C ( fm2/2 , , , i-mij/2\ 
J ai J a2 ' ' ' J J \J a2 ' ' ' J } 



This implies that n — ai is even, and 02^ + • • • + Ofe^ is a partition of (n — ai)/2. Setting 
2i = n — fli , we get 

Jn— /_^Jai Ja2 " ' J ~ J n-2i ^ J ai ' " J Uk ^ J ai ' " J ■ 

7r(n) 0<2<n/4 7r(i) 7r(n) 

ai<n/2 

Using Lemma Ml ) again, this becomes 

fn= ^ /n-2i/i + G(/i, . . . , /l„/2j) 
0<i<n/4 

for a specific function G. ■ 



3 Random power series 

In this section, we consider the reciprocal of a random power series in F2[[q]]. The results of this 
section are strong evidence that the density of J- plays little to no role in determining the density 
of and that unless the coefficients of J- have some structure, the density of T is 1/2. 
Recall that a Bernoulli variable is a random variable that takes values in {0, 1}. 

Theorem 3.1. Suppose that /i,/2, ■ • • are independent Bernoulli variables, with 
infmin{P [/„ = 0] , P [/„ = 1]} > 0. 

n 

Then 6{J-) = 1/2 with probability 1. 



We need the following two lemmas. 
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Lemma 3.2 (Levy's Borel-Cantelli lemma). Let Ei,E2,---, be events, and define Zn '■= 
Sfe=i ^Bfc ' random variable that records the number of Ei, E2, . . . , En that occur. Define 

:=P[Sfc I Ei,E2,...,Ek-i]. 

V'^'kLi^k diverges, then Zn is asymptotically equal to X]I-=i ^fc with probability 1. 

For a proof of Levy's Borel-Cantelli lemma, we refer the reader to 14, Sec 12.15]. 

Lemma 3.3 (Binary Central Limit Theorem). Let Xi be with probability 7^ and 1 with 
probability 1 — 7^, and suppose that Xi,X2, ■ ■ ■ are independent. Then, as n ^ 00, 



^ = (mod 2) 



if and only if some 7^ = 1/2 or X]r=i n^i^lTi' 1 ^ 7*} diverges. 

Proof. Let Sn ■= ^^^"^ define pi by P [Sn = (mod 2)] — pi. Clearly Sn is even if and only 

if Sn-i and Xn arc both even or both odd: 

Vn = Pn-lln + (1 -p„-l)(l - 7„). 

Clearly 2pi — 1 = 271 — 1, and 

2p„ - 1 = 2 (p„_i7„ + (1 -p„_i)(l - 7„)) - 1 = (2p„_i - 1)(27„ - 1), 
which provides the base case and inductive step for the equality 



2pn-l = n(27» - 1) 



i=l 



By the standard results for infinite products, we now see that 2p„ — 1 if and only if 27.; — 1 = 
for some i or X]r=i niii^{7ij 1 ~ 7j} diverges. ■ 

Proof of Theorem \H.l\ We begin with some notation: 
an := P [/„ - 0] , 



min{Q!„, 1 - an}, 
inf /3„, 



and note that < /3 < 1/2. We will show first that P [/„ = O] ^ 1/2 as n 00, and then will show 
that (5(F) = 1/2 with probabihty 1. 



LemmalOIi) says that /„ = /„ + fjfn-j, whence 



[/n=0] =P[/„=0]] 



fjfn-j — 1 



is a weighted average of P [/„ = 0] = a„ and P [/„ = 1] = 1 — a„. Consequently, P [/„ = O] > Pn > P 
and P [/„ = 1] > Pn > P- 
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Set Bn :— {i: < i < n/4, fi = 1}, and set Gn :~ G{fi, . . . , /l„/2j), where G is the function 
from Lemma l2.2f to). We have, from Lemma l2.2r if 

fji — fn-2i + Gn- 

From the previous paragraph, we know that E is at least J^lZo^^ Pi — 1*^/4] • In particular, 

a routine calculation shows that \Bn\ — > oo with probability 1. Thus, P [|i3„| > Kn] ^ 1 if Kn goes 
to infinity sufficiently slowly. We have 



P [/„ = 0] = 

P [fn = I \Bn\ > Kn] P [|S„| > Kn] + P [/„ = I \Bn\ < Kn] P [|B„| < Kn] , 

which for large n becomes P [/„ = O] = P [/„ | |i3„| > Kn] ■ 

We now observe that /„ = if and only if Gn — X^ies /n-2i (call this sum (t„), so that 

F[fn=0\\Bn\> Kn] = 

P [G„ = I |B„| > Kn,an = 0] P [an = | |S„| > Kn] + 

P [Gn = 1 I \Bn\ > Kn, f7„ = O] P [a„ = I | |S„| > Kn] . 

This is a weighted average of P [an — | \Bn\ > Kn] and P [an — 1 | \Bn\ > Kn] , both of which go 
to 1/2 as 71 ^ oo by the Binary Central Limit Theorem. Thus, 

P [/„ = 0] « P [/„ = I \Bn\ > Kn] ~ i 

with each of the becoming as n ^ oo. 

Now that we have shown that P [/„ = O] ^ 1/2, we know that E [(5(f',rt)] 1/2, but this 
does not imply that S{F, n) — > 1/2 ever, much less with probability 1. This last step again requires 
the at-least-weak independence of /„ from /i, . . . , /„_i, and the technicalities are handled for us by 
Levy's Borel-Cantelli lemma. 

Let Ek be the event {fk = 0}, and set := P [/^ = | /i, /2, . . . , fn] ■ By the comment above, 
< /3 < ^fc, so X^fc^i — oo. Thus, by Lemma 

with probability 1. For every e > there is an uq such that for all n > no 



n '■ — ' 71 

fc=l fc=l 



These upper and lower bounds on S{F,n) are non-random, so we may take expectations (for large 
n) to get 



1 " 1 1 " 

fe=i fe=i 
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where we have used the hnearity of expectation and the previously proved E [/„ = O] = P [/„ = O] — * 
1/2. This impHes that i ^fe ^ 1/2 also. Consequently, |(SJ| now implies that 

S{F) := lim S{F,n) = l 
with probability 1. ■ 



4 Polynomials 

In this section, we study the reciprocals of polynomials in F2[(7]. The coefficients of such a reciprocal 
are periodic (see Proposition below) , and so the reciprocal has rational density. We also give 
some indication of how the densities of reciprocals of polynomials are distributed, beginning in 
Subsection 14. II In Subsection l4.2l we use the theory of de Bruijn cycles to exhibit an infinite family 
of polynomials whose reciprocals have densities strictly larger than 1/2; in Subsection l4.3l we show 
that if two polynomials have product 1 + (-D > 4), then at least one of them has a reciprocal 
with density at most 1/2. In Subsection 14.41 we show that the reciprocal of an eventually periodic 
sef^ containing is an eventually periodic set containing 0. 

Let J2iZo ^^2* be the binary expansion of n; we define the polynomial Vniq) ■— J2i^o ^i't ^ ^2[<i\- 
Clearly this indexes all polynomials, and the invertible polynomials are precisely those with n odd. 
For a polynomial Q S F2[g], we let Q be the same polynomial with coefficients (all or 1) in Z. For 
instance, Q — Vq^^-^ for every polynomial Q. 

We denote by i{P) the length of the polynomial, i.e., ^{V) — V{1), and by deg(7') the degree of 
the polynomial. Also, ord('P) is the least positive D such that V divides 1 + . It is not immediately 
obvious that ord('P) is well defined: it is for invertible V, and this is the content of Proposition 14.11 
below. For each polynomial V, we define T'* by PV* = 1 + q°^''^^ . We shall see that the properties 
of P and V* axe intimately related ("Propositions 14.21 and I4.5|) . 

If ord(7-') = 2'^''s('P) _ then V is called primitive, and F2[(7]/ (J^) is isomorphic to F2dog(p) , with 
multiplicative generator q. All primitive polynomials are irreducible, but not vice versa; for example 
I + q"^ + and 1 + q + q"^ + q'^ + q"^ are irreducible but not primitive. 

Figure O tabulates properties of Pn for odd n < 256, including factorizations, P*, and densities 
of reciprocals. 

In Figure^ we plot the points {n,6{Pn)) for odd n less than 2^"^. We note that 6{Pn) tends 
to be near 1/2, but is biased toward being below 1/2. This is also suggested, but not proven, by 
Proposition 14 . 51 below . In Proposition 14.41 we give an algebraically-described infinite set of n such 
that d{Pn) > 1/2. Note that S{Pn{q'')) = \5{Pn{.q)), i.e., if Pn is a polynomial in q^, q^, etc, then 
its density is a priori less than 1/2, 1/3, etc. These points have been plotted with squares. 

In Figure El we plot the empirical distribution function of 5{Pn)- The large discontinuities 
near 1/2 mean that these densities occur with large frequency (fully 421 of the 2048 polynomials 
Pi,P3, . . . ,7^4095 have reciprocals with density exactly 1/2). Again visible in Figure is the pref- 
erence of P to have reciprocal with density less than 1/2. The most interesting issue raised in this 
section, which remains unanswered, is to describe the set 

{S{P) : P is a polynomial}. 

For example, is there an n with S{Pn) =3/4? 

'^More precisely, a set whose indicator function is eventually periodic. 
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Figure 3: Properties of P„ for odd n < 256. The "Factors" column records whether Vn is reducible, 
irreducible or primitive. If Vn is reducible, then the factors column evaluates the factors of Vn at 
q = 2. For example, the factors of P245 are given as 3^ • 7 • 11, whence P245 = P^PrPn- 
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Figure 4: The points {n,5{Vn)) with n odd, except (1,0) and (3, 1) 



4.1 Order and Density 

Our first proposition demonstrates that ord(7') is well-defined, and our next proposition shows the 
connection between 5{V), ord(7'), and V* . 

Proposition 4.1. IfV is a polynomial, then ord('P) is finite.'' 

Proof. Let Vi, . . . , Vfc be the irreducible factors of V, and let di be the multiplicative order of q in the 
field ¥2[q]/{Vi). In particular, 1 + q^'^' is a multiple of Vi for each a; G N. Set L :— Icmjdi, . . . , dk} 
and define V* by ViV^ = 1 + q^ and for 1< i < fc by V^V* = l + q 



. Now 



V ■ nK=(i+'?'^)(i+'7'^)(i+9''^)---(i+'z''"'^) 

i=l 

= (1 + q2^)(1 + g2i)(l + q^^) ■■■{! + q^"'"^ 
by repeated use of the children's binomial theorem. 



) = {l + q' % 



We emphasize that, given V and the equality VJ- = 1 + q^ for some D is not uniquely 
determined. For example, V{q)'P{q)J-{q'^) = 1 + q^^. Nor does the proof given above always provide 
the minimal D. 

Proposition 4.2. S{r) = e{P*)/oTd{P). 



Proof. Since V — 



V 



1, we see that the reciprocal of V is periodic with period ordT' (although 



this may not be the minimal period), and in each period has density £('P*)/ ordT-*. ■ 
4.2 de Bruijn cycle algebra 

Our next proposition shows that the reciprocal of a polynomial is a special case of a linear-shift 
register. Fortunately, there is an enormous literature on linear-shift registers (see for example). 

Actually, the proof can be refined to show that ord('P) | 2'^°^^' — 1 if P is irreducible, and otherwise ord(P) = 
2* lcm{ord{Vi), . . . , ord{Vfc)} for some 1 < 2' < fc, where V = Vi--Vk- 
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0.2 0.4 0.6 0.8 1"^ 

Figure 5: The distribution function A(.t) := 2"" ■ #{n: 1 < n < 2^'^,n odd,(5(P„) < x} 

Proposition 4.3. If J- is a polynomial with degree d, then (letting fj — for negative j) 
d 

fn — fjfn-j- 

i=i 

Alternatively, f„ is the constant term of q^" mod J-. 



(6) 



Proof. Since fj = for all j > d, the recurrence © follows immediately from Lemma l^^j). 

Let M be the matrix whose k^^ row is the elementary vector supported in coordinate fc + 1 , for 
fc = 1, . . . , d — 1, and whose last row is the vector (/o, . . . , fd-i), i-e., M is the companion matrix of 
J-. Write c„ for the constant coefficient of q" mod !F. We claim that 



M 



( \ 






\cs+d-l) 




K^s+d) 



(7) 



To see this, let denote scalar projection of elements of ¥2[q\/{J-) onto ^ and let X denote 
multiphcation by q in F2[g]/(J^). Both of these maps are linear, and it is easy to see that Yfe = 
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Yk-iX-^ + fkYo, for 1 < fc < d - 1. Therefore, 

Yo - Yd-iX-^ 

= Yd-2X-^ + fd-iYoX-^ 

= Yd-^X-'^ + fd-2YoX-^ + fd-iYoX-' 

d-l 

= Y,f,Y,X''-^. 

0=0 

Applying Yq to yields Cs+d — Yl'j=o fj'^s+j, which implies 0. Set a„ :— c_„ (define c on 
negative subscripts by using the recurrence). Thus, the sequences (a„) and (/„) satisfy the same 
recurrence, with initial conditions ao = /o = 1, c_i = Ui = fi ~ (for —d < i < 0). ■ 

Our next proposition computes the density of the reciprocal of every primitive polynomial, and 
thereby produces an infinite family of polynomials whose reciprocals have density greater than 1/2. 

Proposition 4.4. If V is a primitive polynomial with degree d, then S{V) — 

Proof. A de Bruijn cycle of order d is a binary sequence {S{n)}'^^^ in which every binary d-word 
appears in a "window" {S{n + 1), . . . , S{n + d)) for some j (indices taken modulo q'^). A reduced 
de Bruijn cycle is a string of length q"^ — 1 which achieves every d-word in some window, except for 
the word 0''. Note that a reduced de Bruijn cycle may always be turned into an ordinary de Bruijn 
cycle by inserting an extra "0" into its longest run of O's. 

If P is primitive, then q is a generator of and it is a classical result that the sequence 
of constant coefficients of the powers of a multiplicative generator yield a reduced de Bruijn cycle. 
Thus, by Proposition 14.31 the first 2*^ — 1 coefficients of V are a reduced binary de Bruijn cycle of 
order d. The reader wishing to explore de Bruijn cycles further can find the basics in [SJHI. 

Since every string except 0'' appears in P, there are exactly 2'^~^ ones in any period. ■ 

4.3 Polynomials with non-high density reciprocals 

We see in Figure0]that polynomials typically have reciprocals with density near 1/2. In FigureEl it 
is apparent that there is a connection between the density of P and P* . Our next theorem elucidates 
the connection. 

Proposition 4.5. If OTd{P) > 4, then mm{S{P), d{P*)} < 1/2. 

This proposition is best possible in that P51P15 — 1 + q^, and S{P5i) — 5{Pi^) = 1/2. 
Proof. Set D ord(7'). We assume without loss of generality that deg(P) < D/2 < deg{P*). If 
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Figure 6: The points {S{Vn),S{P*)) for odd n < 2^2 



deg('P) < 3, then we appeal to the following table of calculations: 



V 




•p* 


5{V*) 


1 


l + q" 


Y^oo „nD 


l/D 


1 + 9 






2/D 






Er=o(9'"^ + 9'"^+') 


l/D 


l + q + q^ 






3/D 



In the case V = l + q + q^ = -f^, we see also that D = (mod 3), and by hypothesis D > 4, so 
that D/3< 1/2. 

We assume now that deg('P) > 3. Since 

P* = =V + q^V + q^^V + q^^V + ■■■ 

1 + q^ 

and deg(7'_) < D, we have d{P*) = i{V)/D. If P has any zero coefficients, then i{V) < deg(P) <D/2 
and so 6{P*) < 1/2. If V has no zero coefficients, then V ={1 + g<i<'s(P)+i)/(i + in ^hich case 
P = (1 + q)/{\ + gdeg(p)+i)^ ggj.jgg i^^g density 2/(deg(P) + 1). Since deg(P) > 3, this 

quantity is < 1/2. ■ 

Corollary 4.6. IfV ^ {1, 1 + g, 1 + g + q^} is a polynomial and V* is primitive, then S{P) < 1/2. 



4.4 Eventually periodic sets 

An eventually periodic set is one whose generating function has the form £{q) + y^^, for some 
polynomials with deg('P) < D, and exactly one of £,V has constant term 1. The finite sets 
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containing are examples. Another example is given by the set N\{n: n = 2 (mod 4)} (which 
has density 3/4), whose reciprocal is the set {1} U {n G N: n is congruent to 0, 2, 5, or 6 modulo 7} 
(which has density 4/7). 

Proposition 4.7. The reciprocal of an eventually periodic set is an eventually periodic set. 

This proposition is essentially the same as that which asserts that rational numbers have even- 
tually periodic decimal expansions. 

Proof. Obviously, the reciprocal of a ratio of polynomials (each with constant term 1) is a ratio of 
polynomials (each with constant term 1). All that we need to observe is that such a ratio TZ/S can 
be written in the form 



S ' 1 + g^' 

with deg(Q) < D. By long division, we can write TZ/S in the form £ + V/S with dcg{V) < deg{S). 

vs' 



But this is the same as £ + -g^, where SS* = 1 + g^, and deg{rS*) < deg{SS*) = D. 



5 The powers of two 

We saw in Section^ that the reciprocal of a polynomial (other than Vi) has positive density. One 
might wonder if the reciprocal of any set with zero density has positive density. Our next theorem 
shows that this is not the case. 

We note the m = 1 case of Theorem 15.11 the reciprocal of Ai = {0} U {2": n G N} is Ai = 
{2" — 1 : n G N}. This is easily proved directly by considering the following sum-preserving involution 
onAixAi. For s, i G N and distinct, define Ai(0, 0) = (0, 0), a*(0, 2*+i-l) = (2*, 2*-l), /i(2^ 2*-l) = 
(2*,2'' — 1), /i(2*,2* — 1) = (0,2*+^ — 1). The existence of this sum-preserving fixed-point-free 
involution proves that every positive integer n can be written in the form a -I- a, where a G Ai and 
a G Ai, in an even number of ways. A similar proof can be given for m — 2, and presumably for 
any m, but quickly grows tedious. We now give an algebraic proof that docs not depend on m. 

Theorem 5.1. Let m > 1. The reciprocal of the set A,,,, := {0} U {2'"": n G N} is the set 

m—l N 

-1+ ^ a;,2^+""': G {0, 1}, f ^ 0, G N I . 

i=0 ■' 

In particular, both 5 (Am, n) and 5 {Arm n) are Om( ^°^" ) ■ 

Proof. Set T{q) — J2n>o l"^"^'"^ ■ the children's binomial theorem J-{q^) — J-iq)'^, and conse- 
quently by induction we see that T{q^"') = T{q)'^"' . 

Now, by the definition of J-{q^ ) = -^(g) + <Z ^md so 

q^Tiqf"' +T{q) 

= (1 + Hi)) [Hi) + Hi? + Hi? + ■■■ + Hq?^-') 

(m-1 \ 
1+ n {^+Hif)\ 

(m-1 \ 
1+ n 
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The series 1+J^(g) is the generating function of {0}U{2"": n e N}, and 1 + 11110^(1 + -^(9^*)) is the 
generating function of { XiS)^ Xi2*+™"' : Xi € {0, 1}, x ^ 6,ni € N}, so this identity is equivalent to 
the theorem. ■ 

The reader may be interested to note that the reciprocal of the extremely thick set N\ {2" : n e 
N} is the thin set {0, 3} U {2" - 1, 2" - 3 : n > 3}, whereas the reciprocal of N \ {4" : n e N} appears 
to have density 1/2. 

Our next theorem shows that the examples given by Theorem 15 . II are extremal. It is impossible 
for a set and its reciprocal to both grow sub-logarithmically. This result was suggested to us by 
Ernest Croot [personal communication]. 

Theorem 5.2. Let F be reciprocals (not both {0} ), and suppose that r is the least positive integer 
inFUF. Then 

\Fn[0,n] \ + |Fn [0,n]| > 2 + [log2(n/r)J. 

Proof. First, note that r e F n F. Let iV > r, so that neither F n [1, N) nor F n [1, TV) is empty, 
and let m, m be the largest elements of those sets. Since 5™+™ occurs in the product J^f" at least 
once, it must occur at least twice. Since N < m + rh < 2N, we see that 

|i^n [A^,2A^)| + |Fn [iV,2iV)| > 1. 

Straightforward counting concludes the proof, since FU F contains twice, and must intersect each 
of the intervals [r, 2r), [2r, 2'^r), [2'^r, 2^r), .... ■ 



6 Theta functions 

Every quadratic that takes integers to integers can be written in the form cq + cin + C2 "("^"-"-^ with 
Ci £ Z. We wish to study the ranges of such quadratics, but we only wish to consider sets that 
contain 0; without loss of generality we may take cq = 0. Thus, we set 

j tiItl — 1 ) I 
6(ci,C2) := < cin + C2 : n € Z > . 

Moreover, we are only interested in those sets that consist of nonnegative integers, so we may assume 
that C2 > Ci > 0. And since 0(ci,C2) = 0(c2 — ci,C2) we may also assume that C2 > 2ci. Finally, 
we are only interested in those sets whose gcd is 1: we can assume that gcd(ci,C2) = 1. The only 
set with ci = not excluded is 0(0, 1) — {(2) : n > 1}, and the only set with C2 = 2ci that is not 
excluded is 0(1, 2) = {n^ : n>0}. Otherwise, we have C2 > 2ci > 0, and gcd(ci, C2) = 1. 

In Figured we give the number of elements in the reciprocal of 0(ci,C2) (with C2 < 18) that 
are at most 10^. We note that none of the entries of this table are larger than 50450, and the entries 
that are less than 49750 are exactly those with C2 = 2 (mod 4). This computation partially justifies 
Coniecture ll.il 

There is another property of 0(ci,C2) that happens exactly when C2 = 2 (mod 4): the set 
0(ciiC2) is not imiformly distributed modulo 4. 

Proposition 6.1. Let gcd(ci, C2) = 1. The set 0(ci, C2) is uniformly distributed modulo every power 
of 2 if and only if C2 ^ 2 (mod 4) . 
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Figure 7: The number of elements < 100000 in the reciprocal of G(ci, C2) 

Proof. First, suppose that ci = 2/e + 1 and C2 = 4£ + 2. Set 

/(n) := cm + 02 "^^"'^^ = {2e + l)n2 + 2{k - t)n. 

If fc and (. have the same parity, then /(n) = (2£+ l)n^ (mod 4), and since in? takes only two values 
modulo 4, the set is not uniformly distributed modulo 4. If k and I have different parity, then 

(2f + l)n^ + 2{k - e)n = {2£ + l)n^ + 2n (mod 4) 

only takes on the values 0,3 modulo 4. Thus, if C2 = 2 (mod 4), then 6(ci,C2) is not uniformly 
distributed modulo 4. 

Now suppose that C2 = 4£, and since gcd(ci, C2) = 1, we know that ci is odd. We have 

f{n) := cm + C2 "^'^ ~ = 2£n2 + (ci - 2e)n = n (mod 2). 

The formal derivative of /(n) is iin + ci —2i ^ (mod 2). By Hensel's Lemma^, the range of the 
polynomial /(n) hits every congruence class modulo every power of 2. Since for every j, f{n) is 
periodic modulo 2^ with period 2^, we see that it is uniformly distributed modulo 2-' . 
Now suppose that C2 = 2£ + 1 is odd. Set 

G := {{2e + l)m(2m - 1) + Ci(2m) : m G Z} 
H := {{2e + l)(2m + l)m + Ci(2m + 1) : m e Z} 

^Hensel's Lemma: If /(n) is a polynomial with integer coefficients, and the two congruences f{n) = a 
(mod p),f'{n) ^ (mod p) have a simultaneous solution, then /(n) = a has a unique solution modulo every power 
of the prime p. 
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so that e(ci,2£ + 1) ^ GU H. The set G is the range of g{m) /(2m) = (2^ + l)m(2TO - 
1) + 2cim = m (mod 2), which has derivative g'(m) = 1 (mod 2), and the set H is the range 
of h(m) := {2£ + l)(2m + l)m + ci(2m + 1) = m + ci (mod 2), which has derivative h'(m) = 1 
(mod 2). Thus, by Hensel's Lemma, both G and H exhaust every congruence class modulo 2^, and 
by periodicity of g{m) and h{m) are therefore uniformly distributed modulo 2^ . ■ 

6.1 The squares 

Let S{q) — X^S^o ' ^"^"i ^ — {Oj 1) 4, 9, 16, 25, . . . }. Figure El shows 8{S, x) for two ranges of x. 
On the small scale, we see that the relative density behaves irregularly, with many small increases 
and decreases. On the larger scale, we see that the relative density seems to decrease inexorably. 

We characterize completely the values of S in the residue classes 0, 1, 2 (mod 4). 

Let i'p{n) be the integer such that p''p(") | n and p''p(")+i \ n, so that 

p prime 

for every n. Let r2(n) be the number of representations of n in the form + z^, where y and z are 
integers. 

Theorem 6.2. Let n G N. If n is even, then n € S if and only if n is twice a square. If n = 1 
(mod 4) is not a square, then n € S if and only if Vp{n) is even for every prime p except one, and 
that prime p and i/p(n) are both congruent to 1 modulo 4- If n = \ (mod 4) is a square, then n £ S 
if and only if Vp{n) = 2 (mod 4) for an even number of primes p = 1 (mod 4). 

We will need the following lemmas. The first expresses s„ in terms of the number of represen- 
tations of n by a particular (depending on n) quadratic form. The second is quoted without proof 
from '12', and gives a formula for r2{n). 

Lemma 6.3. Let n G N, and let j G N satisfy n = 2^ — 1 (mod 2^+-'^). Then Sn — ^ if and only if 

#|(fco, . . . , : heN,n^ 2^+'fcJ+i + E ^'^^'j 

is odd. 

Proof. By Lemma \2.2( ii). ~ I exactly if there are an odd number of tuples (fco,fci,...) with 
weight 

n = k^ + 2kl + 4:kl + 8kl + -- - . (8) 

Let w(n) be the number of such tuples. We give a weight-preserving involution ji of such tuples, 
and deduce the lemma from 

w{n) = 7^ (fixed points of /i with weight n) (mod 2). 

Since n ^ 2* — 1 (mod 2*+^) for < i < j, reducing © modulo 2,4,..., 2^ successively tells us 
that ko,ki, . . . , kj-i are odd, while n = 2-' — 1 (mod 2-'+^) tells us that kj is even. Now define J to 
be the least integer with the two properties: J > j + 2; and 2k j ^ kj. 

We define 



//(fco, fci, fc2, . . . ) = (fco, ki,. . ., kj-i,2kj, kj+i,k,j, kj,. . ., kj, kj/2, kj+i, kj+2, ■ ■ 
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where kj is repeated J — j — 2 times. That this is a weight-preserving involution is a routine 
calculation. 

The fixed points of ^ are those tuples with = fcj = ~ ^j+a = ■ ■ • • In other words, there 
is a fixed point for each solution to 

71 = fc2 + 2fc2 + . . . + 2-''-ifc|_i + 2^'+ifc|+i. ■ 

Lemma 6.4 ([12, Theorem 3.22]). If Vpin) is odd for any prime p congruent to 3 (modulo 4), 
then r2{n) — 0. Otherwise, r2{n) = '^Y\p{vp{n) + 1), where the product extends over all primes 
congruent to 1 (modulo 4)- 

Proof of Theorem \6.Sl If n is even, then n = 2'^ — 1 (mod 2*'"'"^), so we can apply Lemma l6 . 31 with 
J = to arrive at s„ = 1 if and only if n has an odd number of representations of the form 2kf (with 
ki > 0). Clearly there cannot be more than one such representation, and there is one exactly if n is 
twice a perfect square. 

If n = 1 (mod 4), then we may apply Lemma It) . 31 with j = 1 to arrive at s„ = 1 if and only if 
n has an odd number of representations of the form + 4fc| (with fco and fc2 nonnegative) . 

We assume for now that n is not a square. Since n is odd, there are no such representations 
with fco = 0; E^nd since n is not a square, there are no such representations with fc2 = 0. Thus, every 
such representation fcp + 4fc| gives rise to 8 representations {(±fco)^ + (±2fc2)^, (±2fc2)^ + (±fco)^} of 
n in the form + z^. Moreover, any solution to n — + must have one of j/ or z even and the 
other odd since n is odd, and y =/= z since n is odd. Since n is not a square, neither y nor z is zero. 
Every representation (y, z) occurs as one of a family of 8 such representations, and one of these has 
n = y^ + z^ — y"^ + 4(z/2)^ with y > and z > 0. Thus, s„ = 1 if and only if r2{n)/8 is odd. 

By Lemma 16.41 r2(n)/8 = if Vp{n) is odd for any prime p congruent to 3 modulo 4. Otherwise, 
r2(n)/8 — \ Y\p{vp{n) + 1), where the product extends over those primes that are congruent to 1 
modulo 4 (in the remainder of this paragraph, p is always 1 modulo 4). First, note that i'p{n) is odd 
for some prime p since n is not a square. If some Vp{n) is 3 modulo 4 for some p, then r2(n)/8 is 
even, and similarly if Vp{n) is 1 modulo 4 for two primes p. Thus, r2{n)/& is odd precisely if i'p{n) 
is odd for exactly one prime, and both that prime and Vp{n) are 1 modulo 4. 

Now we assume that n = 1 (mod 4) is a square, say n = . Then, as above, most represen- 
tations of n in the form fcp + 4fc| correspond to 8 representations of n in the form y^ + z'^, but the 
representation n = + 4 • 0^ only corresponds to 4 representations in the form y^ + z'^. Since n is 
a square, we know that Vp^n) is even for every prime p. Thus, s„ = 1 if and only if 

^2(n) -4 
8 

is odd. Using the formula from Lemma [6.41 this happens exactly if 1 = np('^p("-) + 1) (niod 4), 
where the product extends over primes that are 1 modulo 4. This, in turn, happens exactly when 
i>p{n) = 2 (mod 4) for an even number of primes p = 1 (mod 4). ■ 

We suspect that S{S) — and that this may follow from the theory of modular forms, but again, 
this is outside the scope of this paper. We emphasize in CoroUarv 16 . 51 that our characterization of 
S is consistent with Coniecture ll.il 

Corollary 6.5. The set {n Cz N: n ^ S ,n ^ 3 (mod 4)} has zero density. 

Proof. By Theorem 16.21 the set S clearly has no density in mod 2. We will use the description 
given in Theorem 16. 21 to show that S also has zero density in 1 mod 4. 
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By the Wiener-Ikehara Theorem (see W, Section 7.2]), we have for any set A of positive integers 
hm S{A,n) = hm (s — 1) a^'^. 



1 + 

aeA 



Set A — {wp: 1 < n e N,p prime}, and observe that ^(^4) = since 



6iA) < ^hm (s - 1) ^ a-^ 

aeA 



\p prime 

= C(2) (5(primes) = 

Note that the subset of S whose elements are even has density 0, and the subset whose elements 
are congruent to 1 modulo 4 is (except for some squares) contained in A. Thus 

(5({n eS:n^3 (mod 4)}) < (5(squares) + S{A) = 0. ■ 

7 Prouhet-Thue-Morse numbers 

Set t„ = 1 if the binary expansion of n contains an even number of "l"s, and set i„ = otherwise. 
The set T := {n: tn = 1} = {0,3,5,6,9,...} is called the Prouhet-Thue-Morse sequence. This 
sequence frequently arises because it simultaneously has enough structure to analyze, and enough 
"random-like" behavior to be interesting. The survey details four of the occasions that the 
sequence has been independently rediscovered: first in number theory (Prouhet), then combinatorics 
(Thue), then in differential geometry (Morse), and finally chess grandmaster Max Euwe rediscovered 
it to demonstrate that the rules then in use did not imply that chess is a finite game. 

For every n e N, 2n e T if and only if 2n -I- 1 ^ T; thus r{q) E,T=o*"9" ^las 5{T) = 1/2. 
The sequence ioj^ij ■ • • is not eventually periodic (in fact, the real number with binary expansion 
Q.t[)tit2 • ■ • is transcendental [HE]), so T is not a polynomial. A counting argument |2| reveals the 
interesting identity: 

{\ + qfT{qf + {\ + qfr{q)=q. (9) 
Multiplying by f{q) yields gT(g) = [l + q + q^ + q^)'T{q) + l + q^, whence for n > 2 

tn = tn+1 + tn + ^ri-l + 

This leads reasonably directly (albeit with the modest labor involved in deriving 0) to a proof of 
Theorem 17.11 Instead, we give a proof which does not rely on the special form of the functional 
equation lO, and so is more representative of the process of finding reciprocals. 

Theorem 7.1. The reciprocal of the set T of Prouhet-Thue-Morse numbers is 

T — {0} U {4fc ± 1: the binary expansion of k > 1 ends in an even number of "l"s}. 
Consequently, S{T) = 1/3. 
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If (the binary expansion of) k ends in an even number of "l"s, then 4fc + 1 ends with a string 
^Q2fe+i^ (a "1" followed by an odd number of "0"s followed by a single "1"), while 4fc — 1 ends with 
a string 01^*"' (a "0" followed by a positive even number of "l"s). 

Proof. By Lemma l2. II T^q) = X]^o^(^)'?"' where r{n) is the number of ways to write n as 

n = so + 2si + 4s2 + 8s3 H h 2*^5^ H 

where the Sk are Prouhet-Thue-Morse numbers. We will build an involution t on the set of such 
representations, and r{n) will have the same parity as the number of fixed points of t. 

By a tuple, we mean an infinite list of Prouhet-Thue-Morse numbers which is from some point 
on. The weight of a tuple (sq, si, . . .) is X]^o ■^"2". 

We now give the weight-preserving permutation r of the set of tuples which is actually an 
involution. The permutation r has an odd number of fixed points with weight n > if and only 
if the binary expansion of n ends with a string 10^'^"'""'^ 1 (a "1" followed by an odd number of "0"s 
followed by a single "1") or ends with a string 01^*^ (a "0" followed by a positive even number of 
"l"s). These are exactly the numbers of the form 4fc ± 1, where the binary expansion of k ends in 
an even number of "0"s, and this will conclude the proof. 

Defining tiie permutation r: Suppose that sq is even. If sq ^ 2si, then set 

r(so, si, S2, ...):= (2si, so/2, S2, S3, . . . ). 
If So — 2si, then let i be minimal with si 7^ s^, and set 

t(so. Si, S2, ...):= {2Si, Si, Si,. . . , Si, Si, S,;+i, S,;+2, ■ • ■), 

where Si is repeated i — 1 times. The only fixed point with sq even is (0, 0, . . .) with weight 0. 
Now suppose that sq = 3 (mod 4). Since (sq — 3)/2 is even, we can define uq, W2, 1'a, . . . by 

(vq, V2,V3,. . .) := t((so - 3)/2, S2, S3, . . .), 

where the action of r has already been defined above. Note that vi is not defined, and that si has 
not been used. We now set 

t(so. Si, S2, ...):= {2vo + 3, si, t)2, i'3, • • ■ )• 

The only fixed points with sq = 3 (mod 4) are the tuples of the form (3, si, 0, 0, . . . ), where si is a 
Prouhet-Thue-Morse number. These fixed points have weight 3 + 2si. 

Now suppose that sq = 1 (mod 4). If there exists an L such that Si = (sq + l)/2 for < i < L 
and Si = for i > L, then we let t fix the tuple. These will be the only fixed points of r with sq = 1 
(mod 4), and will have weight 2^so + 2^ — 1. Otherwise, if any Si is even (except for the tail of zeros 
in the tuple (sq, si, S2, . . .)), then let K :— min{i : Si even}, and set 

t(so, Si, S2, ...):= (so, si, S2, . . . , sk-i,t{sk, sk+i, sk+2, ••■))■ 

If on the other hand all Si are odd (except for the ending string of zeros), then define vq,vi,V2, ■ ■ ■ 

by 

{vq, Vi,V2, ...):= t(so + 1, Si, S2, S3, . . .) 

and set 

r(so, si, S2, ...):= {vq ~ l,vi,V2, ■ . .). 

That T is an involution with precisely the claimed fixed points is simply a matter of checking 
the various cases; we cheerfully leave this important tedium to the reader. 
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Analysis of r's fixed points with weight n: Suppose that n is even. By parity considerations, 
we see that all tuples (sq, si, . . .) with weight X^i^o ~ have sq even. Since the only fixed point 
with So even is (0, 0, . . .), we see that r(0) — 1 and r{n) is even for all even n > 0. From this point 
on we assume that n is odd. 

Suppose that n = 1 (mod 4), and (sq, si, • • •) is a fixed point of r with weight n. Since n is 
odd, So is either 1 or 3 modulo 4. If sq = 1 (mod 4), then si = 1 (mod 2), and such a tuple can 
be fixed by r only if n = sq, and n is a Prouhet-Thue-Morse number. If sq = 3 (mod 4), then 
si = 1 (mod 2), and such a tuple can be fixed by r only if n = 3 + 2si, i.e., if (n — 3)/2 is a 
Prouhet-Thue-Morse number (and in this case there is exactly one such tuple). Thus r has either 
0, 1, or 2 fixed points, and we care about when it has an odd number of fixed points. Since n = 1 
(mod 4), the binary expansion of n can be written as (xl0'^l)2 for some binary string x and positive 
integer k. We see that the binary expansion of (n — 3)/2 is (x0l'^)2. Thus, if k is even, then either 
both n and {n — 3)/2 are Prouhet-Thue-Morse numbers or neither is. If k is odd, then exactly one 
of n and {n — 3)/2 are Prouhet-Thue-Morse numbers. Hence, t has an odd number of fixed points 
exactly if the binary expansion of n ends in lO'^l, with k an odd number. 

Now suppose that n = 3 (mod 4), and (so,si, . . .) is a fixed point of r with weight n. Since 
n is odd, sq is either 1 or 3 modulo 4. If sq = 1 (mod 4), then si = 1 (mod 2), and such a tuple 
can be fixed by t only if Si — (sq + l)/2 for all < i < L and Si — for i > L. In this case, 
n ~ 2^so + (2^ — 1). Since sq = 1 (mod 4), this implies that the binary expansion of n ends 
with L + 1 "l"s (in particular, at most one value of L can lead to such a fixed point). Moreover, 
2^so + (2^ — 1) is a Prouhet-Thue-Morse number if and only if L is even. If so = 3 (mod 4), then 
si = (mod 2), and such a tuple is fixed if and only if it is of the form (3, si, 0, 0, . . .). This can 
happen exactly if {n — 3)/2 is a Prouhet-Thue-Morse number. 

Suppose that n is a Prouhet-Thue-Morse number. If the binary expansion of n ends in exactly 
2fc > "l"s, then (3, (n — 3)/2, 0, 0, . . .) is the only fixed point of r. If the binary expansion of n 
ends in 2fc + 1 > "l"s, then both (3, {n - 3)/2,0,0, . . .) and 



/n-22'^' + l n+l n+l 




(the term {n + l)/2^'^+^ is repeated 2k times) are fixed points. 

Now suppose that n is not a Prouhet-Thue-Morse number. If the binary expansion of n ends 
in exactly 2fc > "l"s, then 



/n - 22'=-! + 1 n + l n + l 




(the term (n + l)/2^'^ is repeated 2k — 1 times) is the only fixed point. If the binary expansion of n 
ends in 2fc + 1 > "l"s, then there are no fixed points. ■ 
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